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ABSTRACT 

Data from the two MODIS instruments have been accurately geolocated (Earth located) to enable retrieval of global 
geophysical parameters. The authors describe the approach used to geolocate with sub-pixel accuracy over nine 
years of data from MODIS on NASA’s EOS Terra spacecraft and seven years of data from MODIS on the Aqua 
spacecraft. The approach uses a geometric model of the MODIS instruments, accurate navigation (orbit and 
attitude) data and an accurate Earth terrain model to compute the location of each MODIS pixel. The error analysis 
approach automatically matches MODIS imagery with a global set of over 1,000 ground control points from the 
liner-resolution Landsat satellite to measure static biases and trends in the MODIS geometric model parameters. 
Both within orbit and yearly thermally induced cyclic variations in the pointing have been found as well as a general 
long-term trend. 
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INTRODUCTION 

Since February 2000, the first Moderate Resolution Imaging Spectroradiometer (MODIS) instrument [1| has been 
acquiring images of the Earth in 36 spectral bands w'ith nadir spatial resolutions of 0.25 km (2 bands), 0.5 km (5 
bands) and 1.0 km (29 bands). The center spectral wavelengths of these bands range from the visible (0.41 pm) to 
the long-wave infrared (14.2 pm). Two MODIS instruments are currently on-orbit. The first launched on NASA’s 
Earth Observing System (EOS) Terra spacecraft in December 1999 and the second on NASA’s EOS Aqua 
spacecraft in May 2002. Over nine years of data from MODIS/Terra and seven years of MODIS/Aqua data have 
been accurately geolocated to enable accurate retrieval of global geophysical parameters. In this paper, we discuss 
results from the geolocation error analysis approach that has allowed us to maintain and improve the geolocation 
accuracy for both instruments over their mission lifetimes. 

The MODIS geolocation and error analysis approach is discussed in detail in [2] with the first on-orbit results 
presented in [3]. The approach uses a geometric model of the MODIS instruments, accurate navigation (orbit and 
attitude) data and an accurate Earth terrain model to compute the location of each MODIS pixel. The error analysis 
approach automatically matches MODIS imagery with a global set of 1,214 ground control points from the finer- 
resolution Landsat satellite to measure static biases and trends in the MODIS geometric model parameters. Later 
geolocation results are also presented in [4], [5] and [6]. Both within-orbit and yearly thermally induced cyclic 
variations in the pointing have been found as well as a general long-term trend. Also, details of the control point 
matching algorithm are described in [6], 

MODIS data has been reprocessed several times allowing for improvements in instrument characterization and 
algorithms to be incorporated in the standard MODIS products. The current reprocessing started in 2007 and is 
called Collection 5 (C5). The next reprocessing. Collection 6 (C6), is expected to begin in early 2010. For 
geolocation, each collection includes improvements in the geometric characterizations and minor algorithm 
improvements ([4], [5] and [6]). In addition, for C6, an improved Terrain correction algorithm [7] will be added. 
Starting with C5, several improved ancillary' data sets were incorporated: an improved control point library 
containing more than three times the number of control points with a better global distribution [6]; an improved 
digital terrain model based on the Shuttle Radar Terrain Mapping Mission (SRTM) data [8]; and an improved 
MODIS based Land/Water mask. In this paper, we discuss the MODIS/Terra and MODIS/Aqua geolocation results 
from the current C5 reprocessing and the results expected in C6. 



MODIS/Terra Geolocation Results 


The MODIS/Terra C5 geolocation error (Table 1 and Figure 1 ) is small and meets the MODIS geolocation accuracy 
goal of 50 m (1 a), 20% of a MODIS 250 m band. Table I gives the C5 geolocation Root Mean Square Error 
(RMSE) in nadir equivalent units based on an average of 269 control point matchups per day in the track and scan 
(cross-track) direction. Shown in Figure 1 are 


the daily average control point residuals for Table j MODIS control point residuals in nadir equivalent 
each data-day since launch. Also shown are units for C5 and expected C6 results. 
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not follow the prediction. This was corrected 
in the next update that occurred in early 2006. 


An improvement made for MODIS/Terra C5 was the introduction of a sun-angle correction to correct for within- 
orbit variation seen in the ground control point residuals. Because of this correction, as seen in Figure 1, the 
Northern and Southern hemisphere trends are not significantly different from the global trend. 


Based on the C5 results, two improvements will be made in C6. First, a retrospective long-term trend analysis will 
be applied to remove any biases or extrapolation errors. Second, an improved sun-angle correction will be applied 
to correct for residual within-orbit variation. Table 1 and Figure 2 show the expected results from these changes. 
The overall C6 results are expected to improve in the scan direction, reducing the RMSE by 2 m. This is primarily a 
result of the improved retrospective fit in 2005. The long-term trend includes a yearly cyclic term with a peak-to- 
peak amplitude of 26 m in the track and 10 m in the scan direction. This amplitude will be 59% and 29% higher in 
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Figure I . MODIS/Terra C5 control point residuals in the track (top) and scan (bottom) directions. 
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Figure 2. Expected MODIS/Terra C6 control point residuals in the track (top) and scan (bottom) 


the track and scan directions, respectively, than C5. 


The improved sun-angle fit is shown in Figure 3 (black lines). The control point residuals in nadir equivalent units 
are binned in 5° increments (gray filled circles) and displayed with each bin’s standard deviation (gray bars). The 
polynomial fits are used for the sun-angle range from 0° to 120° and then a linear correction for remainder of the 
orbit ( 120° to 360°). The fit has a peak-to-peak amplitude of 40 m in the track and 41 m in the scan direction. This 
amplitude is 43% and 32% smaller in the track and scan directions, respectively, and so compensates for an 
overcorrection that occurred in C5. 
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Figure 3. MODIS/Terra sun-angle fit in the track (left) and scan (right) directions. 


As part of this error analysis, a study was done to detennine whether instrument temperature could be used to 
correct for with-orbit pointing errors since temperature changes are believed to be the primary cause of the 
instrument and spacecraft pointing changes within an orbit. After removing the long-term trend from both the 
residuals and the temperature data from six instrument locations, the temperature that had the highest correlation 
with the control point residuals was the scan mirror temperature (TP SA RCTt MIR thermistor). The relationship 


between this temperature and the pointing error was -5.3 m and -5.6 m per 1° C in the track and scan directions, 
respectively. When this temperature fit was used to remove this wilhin-orbit variation, the overall results were not 
quite as good as the sun-angle fit. Based on this and because of the additional work needed to continuously check 
the quality of the temperature values, the sun-angle fit is the approach that will be used for C6. 

MODIS/Aqua Geoloeation Results 

The MODIS/Aqua C5 geolocation error (Table I and Figure 4) is small and meets the MODIS geoloeation accuracy 
goal in the scan direction and almost meets the goal in the track direction. Table I gives the C5 geoloeation Root 
Mean Square Error (RMSE) in nadir equivalent units based on an average of 232 control point matchups per day in 
the track and scan (cross-track) direction. Overall, the long-term trend is close to zero, except for some periods, 
such as the period in 2007 and 2008, when the overall trend in the scan direction did not follow the prediction. This 
was corrected in the next update that occurred in early 2006. For Aqua, there is a significant with-orbit variation 
that causes the northern and southern hemisphere residuals to vary with respect to one another. In addition, there is 
a small global yearly cyclic term. Neither of these corrections was used in C5 for MODIS/Aqua. 
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Figure 4. MODIS/Aqua C5 control point residuals in the track (top) and scan (bottom) 

directions. 

Based on the C5 results, three improvements will be made in C6. First, a retrospective long-term trend analysis will 
be applied to remove any biases or extrapolation errors. Second, a yearly cyclic fit will be used to remove the global 
yearly variation. The third change will be to implement a sun-angle correction to correct for the within-orbit 
variation. Table 1 and Figure 5 show the expected results from these changes. The overall C6 results are expected 
to improve in both the track and scan directions, reducing the RMSE by 2 m for each. The new long-tenn trend 
yearly cyclic terms will have small, peak-to-peak amplitudes of I m in the track and 5 m in the scan direction. As 
seen in Figure 5, the retrospective fit of the long-term trend that includes the cyclic term bring the global residuals 
close to zero in both the track and scan directions. In (he track direction, almost all of the wilhin-orbit variation is 
removed by the 26 m peak-to-peak amplitude sun-angle correction (Figure 6). There appears to be a small, two- 
cycle per year residual error remaining in the track direction. In the scan direction, the within-orbit variation is 
reduced by the 62 m peak-to-peak sun-angle correction. I lowever, there remains a small north/south difference that 



appears early in each year starting with 2004. Some additional adjustments will also be made in the early post- 
launch period to correct for the initial on-orbit bias. 
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Figure 5. Expected MODIS/Aqua C6 control point residuals in the track (top) and scan 

(bottom) directions. 




The first MODIS/Aqua sun-angle lit is shown in Figure 6. The polynomial fits are used for the sun-angle range 
from 0° to I40 c and then a linear correction for remainder of the orbit (140° to 360°). The MODIS/Aqua scan 
direction sun-angle fit is similar in sign and has a 51% larger amplitude than the MODIS/Terra scan fit. The 
MODIS/Aqua track direction sun-angle fit, however, has the opposite sign and an amplitude that is 35% smaller 
than the MODIS/Terra track fit. 
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Figure 6. MODIS/Aqua sun-angle fit in the track (left) and scan (right) directions. 

As part of the MODIS/Aqua C5 error analysis, a study like the study for MODIS/Terra was done to determine 
whether instrument temperature could be used to correct for within-orbit pointing errors. After removing the long- 


term trend from both the residuals and the temperature data from six instrument locations, there was no temperature 
that had a good correlation with the control point residuals. So, for MODIS/Aqua, no correction based on 
temperature was possible. 


Other Collection 6 Changes 

Several additional changes will be made for geolocation in C6. First, minor improvements in the overall geolocation 
accuracy will be obtained by making small adjustments in several instrument parameters (i.e. yaw bias, pitch bias, 
scan rate). Second, as mentioned above, an improved terrain correction algorithm [7] will be implemented that 
improves 1 km geolocation accuracy in areas with rough terrain. This change will be coupled with a higher quality 
terrain model and an improved land/water mask, both with 0.5 km resolution. Finally, in addition to the 1 km 
geolocation, geolocation will be computed at 0.5 km and stored as offsets from the I km geolocation values. 

Other Geolocation Considerations 

Both the Terra and Aqua spacecraft perform monthly attitude roll maneuvers for calibration purposes. Additionally, 
orbit maneuvers, drag makeup and inclination adjustment, are performed several times a year. For a short time 
before a maneuver and for up to 1.5 hours after a maneuver, the satellite orbit and attitude data are not as accurate, 
so the geolocation has larger errors than mentioned above. Because of this, these maneuver periods are excluded 
when the standard MODIS products are produced in C5. A different approach is being considered for C6 that would 
allow more data near maneuvers to be used. 

For near-real-time users of MODIS/Aqua data, the predicted orbit data provided in the direct broadcast data stream 
is not as accurate as the definitive orbit data that is available 12 to 24 hours later. This is not an issue for 
MODI S/Terra because the direct broadcast data stream contains definitive orbit data. 

Summary and Conclusion 

The error analysis approach described above is being used to accurately geolocate data from both MODIS 
instruments. The current C5 geolocation results were discussed along with the geolocation improvements that are 
expected in the next MODIS reprocessing (C6) set to begin in early 2010. Overall, the accuracy of the MODIS 
geolocation data meets the tough accuracy goal set by the MODIS Science Team. The consistent sub-pixel 
geolocation accuracy enables the accurate retrieval of geophysical parameters globally. Lessons learned from this 
error analysis approach is being used for the Visible Infrared Imager Radiometer Suite (V1IRS) instruments aboard 
the National Polar-orbiting Operational Environmental Satellite System (NPOESS) Preparatory Project (NPP) and 
will be used for future operational NPOESS missions. 
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